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ABSTRACT 



Aims. The chemically peculiar B star if> Phe was, until very re cently, considered a tr iple system, even though the data were not 

conclusive and the orbits rather uncertain. Very recent results by Korhonen et al. (2013) provided a revised orbit, different from the 

then available astrometric Hipparcos orbit. 

Methods. Additional spectroscopic data, obtained with the BESO spectrograph at Cerro Armazones, confirm the newly found orbit, 

even though the resulting radial velocities do not allow to improve on the recent orbit. We combine the latter with the Hipparcos 

measurements to secure the astrometric orbit, and derive the inclination of the system. Using evolutionary tracks, we can finally 

constrain all the parameters of the two components in this system. 

Results. We confirm the mass of the primary, 3 M , and find that the companion has a mass of 0.9 M . The inclination of the system 

is i = 93° + 4.7°, and is potentially eclipsing; we predict the time of the next conjunction. Given that the eccentricity of the orbit and 

the exact value of the semi-amplitude of the radial velocity relies on just one set of points, we also urge observers to measure radial 

velocities at the next periastron passage in April 2015. 
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1. Introduction 

Phe (HIP 8882, HD 11753, HR 558) is a 5th magnitude 
late-type B star whose chemical peculiarity was first noticed 
by A.R. Hyland (Dworetsky 1969). Its MK classification is 
B9pHgMn (Bus combdll984h . 

It was suggested very early on (Campbell & Moore 1928) 
that the radial velocity of (p Ph e was varying. However, 
that claim was only confirmed by Dworetskv et al. (1982) on 
a dataset too small to y ield any orbit, not even the period. 
iLeone & Catanzarol (1199 9) completed the dataset and derived 
the first spec t roscop ic orbit (<? = 0.32, P = 41.489 days). 
iGrenier et al.l (1 1999b also complemented the dataset with one 
more radial velocity but did not derive any orbit. In the original 
reduction of the Hipparcos data (ESA 1997), a circular nearly 
edge-on 878-day orbit was derived from the sole astromet ric ob- 
serva tions. <p Phe thus made it to the triple star class (Tokovinin 
2008) despite the excellent astrometric fit (F2 = 0.10), leaving 
no room for any detectable wobble caused by the inner compo- 
nent. 

Even thou gh the spe ctroscopic orbit has lately been substan- 
tially revised (iKorhonen et al. 2013) with a period now reaching 
sightly more than three years, the multiplicity of the system was 
not questioned. Using the Hipparcos data, we are going to show 
that <p Phoenici turns out to be a genuine binary (Sect. [3}, thus 
allowing the system to be fully characterised (Sect.|4j>. 
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2. Revised spectroscopic orbit 

Using about 150 Coralie spectra taken in 2000, 2001, 2009, and 
2010 together with a few FEROS and HARPS data poin ts on 
top of the older radial velocities, IKorhonen et al.1 (|2013) have 
recently drastically changed the orbit. Instead of a month-long 
period, they obtained a 3-year orbit. Although these new data 
strongly constrain the top of the velocity curve, its bottom (and 
thus partly its semi-amp litude) is still essent i ally se t by the sole 
six additional data from lLeone & Catanzarol (I 1999b . 

The Bochum Echelle Spectroscopic Observer (BESO) is a 
fibre-fed spectrograph based on slightly improved blueprints of 
the European Southern Observatory FEROS instrument (located 
at La Silla). It belongs to the Asttonomical Institute of the Ruhr- 
Universitat Bochum (AIRUB) and is located at Cerro Armazones 
in the Atacama desert in Chile through a partnership between 
AIRUB and the Universidad Catolica del Norte in Antofagasta. 
BESO is attached to the 1 .5-m Hexapod Telescope. It operates in 
the 370-860 nm wavelength range with an average resolution of 
48 00 0, varying from order to order between 43 000 and 60000 
(see lFuhrmann et aDl201 ll for more details about BESO). 

For purely astrometric reasons, <p Phe has been monitored by 
BESO since late 2011. The initial goal of that monitoring was 
to secure the spectroscopic orbit of the inner componen t s (eve n 
the period was still uncertain after Leone & Catanzarol (|T999)) 
which could then be used to improve the Hipparcos solution. 
Seventeen spectra were obtained between November 2011 and 
October 2012 yielding the radial velocities listed in TabQ] The 
radial velocities were determined both by fitting some of the 
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Table 1. BESO radial velocities. Typical uncertainty is 2 km s 



Date 



HJD-2400000 Vel(kmsl) Date 



HJD-2 400 000 Vel (km s" 1 ) 



2011/11/09 


55875.52796296 


6.52 


2012/09/07 


56178.85435184 


16.13 


2012/06/22 


56101.90094908 


10.23 


2012/09/16 


56187.67648149 


13.757 


2012/06/27 


56106.8688426 


10.63 


2012/09/27 


56198.61953704 


9.19 


2012/07/06 


56115.93180557 


7.36 


2012/10/14 


56215.67843749 


12.93 


2012/07/15 


56124.86010416 


10.85 


2012/10/16 


56217.65494213 


13.245 


2012/08/10 


56150.84261572 


12.87 


2012/10/20 


56221.62182871 


17.56 


2012/08/15 


56155.89633103 


10.15 


2012/10/24 


56225.60006944 


15.30 


2012/08/29 


56169.83555555 


13.485 


2012/12/04 


56266.68113426 


15.54 


2012/09/02 


56173.69356481 


13.952 





















strongest lines and by cross-correlation using the xcsao task in 
IRAF. Depending on the quality of the spectrum - the S/N varies 
between 30 to 70 - the formal errors of the cross-correlation is 
between 0.3 and 1 km s . The real error is higher as it combines 
systematics from the wavelength c alibration and intrinsic vari- 
ability of the object as reported by Korhonen et al. (2013). The 
resulting rms for these sole data is 1.92 kms -1 , likely caused 
by the BESO calibration itself which was reported to have had 
some problems during our run. At the time our invest igation on 
</> Phe was initiated, the results from Korhone n et al.l had not 
been released and the adopted S/N was thought to be enough to 
assess the 41 -day orbit. However, unlike all the other existing 
datasets, this one covers 4 consecutive months (as only one mea- 
surement was taken in November ' 1 1 before <p Phe disappeared 
for 6 months). 

Despite the extended phase coverage, the precision of these 
BESO radial veloc ities is not goo d enough to substantially re- 
fine the orbit from Korhonen et al.l (120131) . However, given their 
continuity, our data point s definitively exclude any 4 0-day pe- 
riod orbit as proposed by iLeone & Catanzarol (1 1999b . The or- 
bit from Korh onen et al.l is plotted in the left panel of Fig. Q] 
The scatter of both the BESO data and the old velocities from 
Campbell & Moore (1928) make the usefulness of these sets 
rather limited, especially the latter as far as constraining the pe- 
riod is concerned. 



3. Astrometric counterpart 

The original choice of adopting a circular orbit to fit the Hippar- 
cos data of this object was for the sake of simplification. Indeed, 
it means that not only the eccentricity but also the argument of 
the periastron (ai) can be set to 0. That orbit is listed in the first 
column of Tab. [2] Strangely enough, the very same orbit was 
als o adopted in the new reduction of the Hipparcos observations 
by Ivan Leeuwenl (120071) despite the seemingly poor fit resulting 
from that choice: 0.1 versus 12.34 for the goodness of fit (Wilson 
& Hilferty's cube root tr ansformation which follows a N(Q, 1) 
distribution (IStuart & Ord 1994)). However, in this second re- 
duction of HIP 8882, there was no new orbit fitting per se, the 
original orbit was simply adopted. 

In their general attempt to reprocess the Hipparcos observa- 
tions of the spectroscopic systems listed in S 9 B (Pourbaixet al. 
2004), Jancartetal. (2005) concluded that the orbit by 
Leone & Catanzarol (11999 ) was not present in the astrometric 
data. Neither the spectroscopic orbit nor the astrometric one 
were ever questioned despite the absence of hint of the latter 
in the radial velocities. 

The astrometric orbits resulting from both the Campbell and 
the Thiele-Innes' approach are listed in Tab. [2] The assess- 
ment of such spectro-astrometric combinations was extensively 



described and improved in sev e ral pap ers (Pourbaix & Arenoul 
120011: iPourbaix & Boffinl 120031: Uancart et al.l 120051) . It essen- 
tially consists in fitting the astrometric data with two different 
approaches (freezing some parameters in one while leaving them 
free in the other) and in evaluating the consistency between these 
two solutions through several statistical tests. In both cases, the 
eccentricity, period and periastron time are adopted from the 
spectroscopic orbit. In Campbell's approach, the amplitude of 
the radial velocity curve and the argument of the periastron are 
also adopted from the spectroscopic side. In the alternative, the 
four Thiele-Innes constants are left free. In no case, the radial 
velocities are fitted so the K\ listed in Tab. [2] is derived from 
the astrometric solution alone. If both solutions are consistent 
(which is the case here), Campbell's one should be favoured as 
it results from the larger number of degrees of freedom. 

All the assessment indicators are listed in Tab.|3]and confirm 
the excellent agreement between the spectroscopic orbit and its 
astrometric counterpart. The goodness of fit is significantly re- 
duced with the eccentric orbit with respect to the original circular 
one. Despite the change of the eccentricity and period, the par- 
allax and the proper motion, last three lines of Tab. [2] remain 
unchanged (within the error bars). 

Table 3. Assessment of the presence of the spectroscopic or- 
bit in the Hipparcos ob servations. The notations are those from 
Pourbaix & Bofrfl J2003h 



Indicator 


Value 


Indicator 


Value 


p n 


< 10" 6 


e 


0.589825 


Pri 


< io- 6 


D 


1.922401 


Pn 


< io- 6 


Pr 4 


0.645796 


F2-YI 


-0.866067 


Pr 5 


0.408080 



Although the Thiele-Innes solution fits the Hipparcos data 
very well, the resulting K\ is s lightly larger than the spectro- 
scopic value of lKorhonen et alj d2013l) . i.e. 9.21 kms -1 . An as- 
trometric value smaller than its spectroscopic counterpart is of- 
ten the sign that the assumption about the absence of light from 
the secondary is questionable (very unlikely, as we will see at 
the end of Sect. @}. In the case of HIP 8882, could the spectro- 
scopic value be underestimated? This cannot be ruled out as the 
bottom of the velocity curv e is cons trained by the six points over 
five consecutive days from lLeone & Catanzaro (1999), with the 
observed velocity ranging from -4.4 to -3.0 kms-'. The next 
minimum of the radial velocity will take place on JD 2457 130 
(2015.290). 

For the sake of completeness, we did fit the data com- 
ing from the second reduction of the Hipparcos observations 
(Ivan Leeuwenl l2007). The resulting «o is way too small (4.20 + 
0.33 mas) and, even though the goodness of fit goes from 12.34 
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Table 2. Astrometri c orbits based on the original Hipparcos data with and without spectroscopic help. The second column lists the published 
solution dESAll99 7 f). As in the Hipparcos Catalogue, F2 denotes the Goodness of Fit and m the parallax. The third column lists the proper motion 
after Tvcho-2 ( fog et al.ll2OO0h . 



Element 



Original 



Spectro-based 
Tycho-2 Campbell Thiele-Innes 



e 


(fixed) 






0.589 (fi 


xed) 


P(d) 


877.6541 + 117.3686 






1126.11 (fixed) 


r (JD-2400000) 


48 589.2621 + 38.5598 






53 766.2 (fixed) 


flo (mas) 


5.91+0.93 






8.4 + 0.52 


9.8 


i (deg.) 


93.51+5.34 






93 ± 4.7 


93 


co (deg.) 


(fixed) 






201.6 (fixed) 


196 


Q (deg.) 


161.88 + 6.94 






159 + 6.4 


156 


A (mas) 










8 + 1.0 


B (mas) 










-4± 1.0 


F (mas) 










-2.7 + 0.80 


G (mas) 










-0.6 + 0.67 


F2 


0.10 






-0.83 


-0.87 


K x (kms" 1 ) 


6.93 






9.21 (fixed) 


11 


m (mas) 


10.55 + 0.69 






10.87 + 0.67 


10.6 + 0.69 


/j m (masyr 1 ) 


-34.32 ± 0.68 


-33.7 + 


0.7 


-35.3 + 0.74 


-35.7 + 0.78 


/j s (masyr 1 ) 


-28.34+1.39 


-26.0 + 


0.6 


-25.8 + 0.52 


-25.4 + 0.61 
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Fig. 1. Left panel: Spectroscopi c orbit with data fromlLeone & Catanzarol (1 1 9991) . iGrenier et al J d!999h . iKorhonen et all d2013h (all as filled 
disks), BESO (open triangles) and Campbell & Mooig dl928t) as filled diabolos. Right panel: New astrometric orbit based on the original 
Hipparcos observations. 



to 5.39, all the statistical indicators are made useless by the poor 
weighting scheme adopted in that reduction. It is therefore safer 
not to use it, at least for binaries. 

Although a perfectly edge-on orbit would be consistent with 
the astrometric solution, there is n o sign of an y eclipse ever re- 
ported. The Hipparcos photometry (ES All997l) shows a constant 
brightness during the whole mission. The difference between the 
5th and 95th percentiles is 0.02 Hp mag and there is no sign of 
brightness decrease at the epoch of conjunction. 



Finally, the revision of the orbit leads to a significant change 
in the proper motion with respect to the Hipparcos one, es- 
pecially along t he declination ax is, getting much closer to the 
Tycho-2 value (H0 g et al]|2000h . Unfortunately, the right as- 
cension componen t moves in the opp osite direction. How- 
ever, the UCAC4 dZacharias et al.ll2013i) proper motion (-33.9 + 
1.0; -28.2 ± 1.0)masyr _1 confirms neither the Tycho-2 one nor 
ours, illustrating once more our poorly known this bright system 
is. 
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4. Astrophysical outcome 

Besides the genuine binary nature of <p Phe, our new combined 
astrometric-spectroscopic solution confirms the parallax of the 
system and the orbital inclination, which is now safe to combine 
wit h the mass f uncti on to derive the mass of the secondary. 

iDolk et all d2003l) reported T eff = 10,612 ± 200 K and 
log<7 = 3.79 + 10 tog ether with [Fe/H] = 0.15 from 
Smith & Dworetsk vl (Il993l). On t he grid of stellar evolution- 
ary tracks from Bertelli et al. (2009), the 3-solar mass track with 
Y = 0.3 and Z = 0.017 gives the best match as far as loggr 
and r a ff are concer ned (left panel of Fig. |2). According to 
iHakkila et ail d!997l) . the interstellar extinction A v = 0.03 +0.17 
which coupled to V - 5.112 yield the absolute magnitude 
My — 0.26 ± 0.13. The point (T e g, My) is in ex cellent agree- 
ment with the logage(yr) = 8.41-isochrone from (iBertelli et al.l 
2009) as illustrated in the right panel of Fig. [2] 

From the spectroscopy, f(m) - 0.048 + 0.0015M o , which, 
combined with i = 93° (astrometry) and M\ - 3.0 + O.12M 
(evolutionary track), lead to M 2 - 0.91 + O.O25M . Using 
one solar mas s as initial gue ss for M%, the dynamical paral- 
lax method (lBinnendiiklll960l) yields M 2 = O.9O45M after a 
few iterations, very consistent with our track estimate. On the 
logage(yr) = 8.41-isochrone (Pig. [2j, such a mass corresponds 
to the big point on the right panel. The magnitudes of the two 
components differ by 5.7 in V and 3.9 in K. Such a large differ- 
ence in V explains the absence of signal in the Hipparcos data 
and in the spectra. The semi-major axis of the relative visual or- 
bit is 36.3 mas only, making this sy stem essentially inaccessible 
to VLT/NACO (ISchoUeretal.ll2010t) . The two components are 
nevertheless bright enough to be resolved with the VLTI. 

Even though, with Am - 5.7 in V, some hints of the sec- 
ondary peak might be detectable in FEROS-based cross corre- 
lation function, a thorough inspection of the publicly available 
spectra has not revealed anything. Should eclipses occur (JD 
2457 078, i.e. 2015.148), they would yield a change in magni- 
tude of 0.010 and 0.057 in V and K respectively, so the constant 
brightness reported by Hipparcos cannot be used to rule out the 
possibility of eclipses. Such eclipses would also constrain the 
orbit to be edge-on within 0.0001°. 
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5. Conclusion 

It is rather disappointing to note that despite the fact that cp Phe 
belongs to the Bright Star Catalogue, in which it was noted for 
its chemical peculiarities, one has to wait till 2013 to have the 
first plausible characterisation of the components of this binary 
system. 

Out of the 235 orbital solutions published in Hipparcos (ESA 

1997), 115 were derived without using any ground-based solu- 
tion, 70 of which had their eccentricity and periastron angle ar- 
bitrarily set to as for <p Phe. How many other systems deserve 
some revision? 
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Fig. 2. Left panel: Four evolutionary tracks from dBertelli et a l. 2009). Right panel: logage(yr) = 8.41-isochrone with same Y and Z as in the 
best evolutionary track. 
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